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ABSTRACT 

We present the spectral energy distribution (SED) of the pre-planetary neb- 
ula, IRAS 19475+3119 (1 19475), from the optical to the far-infrared. We identify 
emission features due to crystalline silicates in the ISO SWS spectra of the star. 
We have fitted the SED of 1 19475 using a 1-D radiative transfer code, and find 
that a shell with inner and outer radii of 8.8 x 10^^ and 4.4 x 10^ ''cm, and dust tem- 
peratures ranging from about 94 K to 46 K provide the best fit. The mass of this 
shell is >1 [34 cm^ g-V/t(lOO/im)][(5/2OO]M0, where ^(lOO/im) is the lOO/xm dust 
mass absorption coefficient (per unit dust mass), and 5 is the gas-to-dust ratio. 
In agreement with results from optical imaging and millimeter-wave observations 
of CO emission of 1 19475, our model fits support an r~^ density law for its dust 
shell, with important implications for the interaction process between the fast 
collimated post-AGB winds and the dense AGB envelopes which results in the 
observed shapes of PPNs and PNs. We find that the observed JCMT fiux at 
sub-millimeter wavelengths (850/im) is a factor ~2 larger than our model fiux, 
suggesting the presence of large dust grains in the dust shell of 1 19475 which are 
not accounted for by our adopted standard MRN grain size distribution. 

Subject headings: Stars: Circumstellar Matter, Stars: AGB and post-AGB, Stars: 
Individual: Alphanumeric: IRAS 19475+3119 



1. Introduction 

Pre-Planetary nebulae (PPNs), objects in transition between the AGB and planetary 
nebula (PN) evolutionary phases, hold the key to our understanding of the relationship 
between these late evolutionary stages of low and intermediate mass (~ 1 — 8Mq) stars. 
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The hydrodynamic interaction of one or more fast, coUimated post-AGB outflows, believed to 
occur during the PPN phase, with the dense, massive, slowly-expanding spherical AGB wind 
is believed to be responsible for shaping planetary nebulae (Sahai & Trauger 1998, Sahai 
2002). Imaging surveys with the Hubble Space Telescope have been crucial in revealing a 
wide variety of aspherical morphologies in PPNs (e.g. Sahai 2004; Ueta et al. 2000), however 
the bulk of the circumstellar mass, which often resides in a spherical, dusty component 
surrounding the central aspherical nebula, is not well-probed in these data. Thus a crucial 
ingredient for theoretical studies of such interactions, namely the mass and density of the 
ambient circumstellar medium (e.g., Lee & Sahai 2003, 2004), is lacking for a majority of 
the PPNs discovered in HST surveys. Although attempts have been made to estimate the 
mass from submillimeter data (e.g. Glcdhill ct al. 2002), these are compromised by the 
simplifying assumptions about the dust temperature. 

The mineralogy of circumstellar dust shells around PPNs became possible with the 
wavelength coverage and spectral resolution offered by the Infrared Space Observatory (ISO) 
(see e.g. Hrivnak et al. 2000; Molster ct al. 2002a, b, c; Gauba & Parthasarathy 2004 and 
references therein). In order to understand the grain temperatures, mass-loss history and 
morphology responsible for the formation of these minerals in the circumstellar environment 
of PPNs, detailed modelling of the spectral energy distributions (SEDs) of a few PPNs (e.g., 
HD161796, Hoogzaad et al. 2002; IRAS 16342-3814, Dijkstra et al. 2003a; IRAS 16279-4757, 
Matsuura et al. 2004; Gauba & Parthasarathy 2004) have been undertaken. 

In this paper, we examine the composition and distribution of dust in the PPN, IRAS 19475+3119 
(hereafter 1 19475) and model its SED from the optical to the far-infrared to set constraints 
on the physical properties of its dust shell (mass, temperature, size). 119475 is listed as a 
F3Ia star in the Case-Hamburg luminous stars survey (Stock et al. 1960). High resolution 
optical spectra of the star have confirmed its post-AGB nature (Klochkova et al. 2002; 
Sivarani et al. 2001; Arellano Ferro et al. 2001). A point-symmetric nebula with an extent 
of 4.9" X 3.4" has been seen in J-band polarized flux images (Gledhill et al. 2001). CO J=l-0 
and 2-1 emission typical of AGB and post-AGB objects was detected from the circumstellar 
envelope of the star (Loup ct al. 1993). Hrivnak & Bieging (2005) noticed asymmetry and 
possible structure in the CO J=2-l and 4-3 lines. Recent imaging of 1 19475 with the Hub- 
ble Space Telescope at optical wavelengths reveals a dusty quadrupolar nebula and spherical 
halo, surrounding the central star (Sahai 2004; Sahai et al. 2006). 

The rest of the paper is organized as follows. In §2, we provide observational details 
of the ISO data, in § 3 we describe the ISO spectra, and the re-construction of 1 19475's 
spectral energy distribution (SED) from these as well as broad-band photometric data, in 
§ 4 we describe our modeUing of the SED using a spherical dust radiative transfer code, in 
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§ 5 we discuss our models, and in § 6, we present our conclusions. 

2. ISO Observations 

Infrared Space Observatory (ISO) observations of 1 19475 were extracted from the ISO 
data archive. These include spectroscopic observations made with the Short Wavelength 
Spectrometer (SWS) and the Long Wavelength Spectrometer (LWS) and spectrophotometric 
observations made with the imaging photopolarimeter (ISOPHOT) onboard ISO. An off- 
source measurement made with the LWS in order to estimate the background was also 
extracted from the archive. A log of the observations is given in Table 1. 

ISO SWS spectra have a wavelength coverage of 2.38—45.2 //m. Our spectra were 
obtained in the low resolution mode (AOT 01) of the SWS instrument (de Graauw et al. 
1996) with a 33" x 20" aperture. Each SWS spectrum contains 12 subspectra, that each 
consist of two scans, one in the direction of decreasing wavelength ('up' scan) and the other 
in the direction of increasing wavelength ('down' scan). There are small regions of overlap 
in wavelength between the sub-spectra. Each sub-spectrum is recorded by 12 independent 
detectors. 

LWS spectra extend from 43— 197/im. The LWS observations were obtained in LWSOl 
mode, covering the full spectral range at a resolution (A/AA) of ~ 200. The LWS circular 
field of view had a diameter of 84". The ISO LWS instrument and its calibration are described 
in Clegg et al. (1996) and in Swinyard et al. (1996) respectively. 

ISOPHOT observations were carried out using the spectrophotometer subsystem PHT- 
S. PHT-S consists of two low-resolution grating spectrometers covering the wavelength ranges 
2.5—4.9 iJ,m (PHT-SS) and 5.8—11.6 /^m (PHT-SL) and having a common entrance aperture 
of 24" X 24". Each channel has a hnear 64-element array of Si:Ga detectors. The spec- 
tral band width (FWHM) of a single detector is 0.0383 //m for PHT-SS and 0.0918 //m for 
PHT-SL, resulting in a mean A/AA of about 95 for both channels. A more detailed de- 
scription of PHT-S has been given by Klaas et al. (1997). In the PHT40 observing mode 
spectrophotometry is performed simultaneously at wavelengths 2.5—4.9 /im and 5.8—11.6 
/im. 
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Table 1: Log of ISO observations 



Object 


Instrument 


Date of Obs. 


Duration of Obs. (s) 


TDT^ 


Mode^ 


Speed'^ 


IRAS 19475+3119 


sws 


12/11/1996 


1140 


36100905 


SWSOl 


1 


IRAS 19475+3119 


sws 


19/04/1997 


6538 


52000931 


SWSOl 


4 


IRAS 19475+3119 


LWS 


12/11/1996 


1266 


36100904 


LWSOl 




Off-source 


LWS 


12/11/1996 


1268 


36100943 


LWSOl 




IRAS 19475+3119 


PHT-S 


19/04/1997 


364 


52000932 


PHT40 





^TDT number uniquely identifies each ISO observation. ^Observing mode 
used. ^Speed corresponds to the scan speed of observation. 
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3. Analysis 

Details of data reduction are described in Appendix A. The reduced SWS, LWS, and 
PHT-S spectra are shown in Figs.l, 2 & 3 , respectively. 

3.1. Infrcired spectral features 

Emission features at 33.6 /im and 43 //m were observed in the SWS spectra (Fig. 1). 
The 33.6 /xm feature is attributed to forsterite, a crystalline silicate (Waters et al. 1996; 
Waters & Molster 1999; Gauba & Parthasarathy 2004). Both pyroxenes (such as enstatite) 
and water-ice show features at ~ 43 /xm. ISO observations of post-AGB stars have shown 
that the 43 /im feature is almost always accompanied with a prominent emission feature at 
40.5 iim attributed to pyroxenes such as enstatite and/or diopside (Molster et al. 2002b). 
The non-detection of the 40.5 //m feature in 1 19475 is therefore puzzhng, suggesting that 
pyroxenes are not (or only very weakly) present in this object 

The only other post-AGB object which shows a spectrum in the ~40 //m region similar 
to that of 1 19475, is HD161796 (Molster et al. 2002b), which, like 1 19475, is a high latitude 
F3Ib post-AGB star with a detached circumstellar dust envelope (Parthasarathy & Pottasch 
1986; Hoogzaad et al. 2002). HD161796 shows a prominent 43 iim feature and a weak 
40.5 /xm feature - in order to fit these features, Molster ct al. (2002c) need a combination 
of pyroxenes and crystalline water ice, with the latter dominating the ~40 emission 
complex. This is because in the ~40 fim region, crystalline water ice contributes a feature 
only at 43 fim, hence its dominant presence helps to reduce the intensity of the 40.5 /im 
feature (produced by pyroxene only) relative to the 43 /xm one (produced by pyroxene and 
ice). The presence of crystalline water ice in HD161796 is supported by the prominent 62 
/im feature in its LWS spectrum. 

Similarly, if the 43 /xm feature in 1 19475 is due to (or dominated by) crystalline water 
ice, then one would expect to see a crystalline water-ice feature at ~ 62 /im. We now examine 
the possibility that the latter is not seen because it may be much weaker compared to the 
43 /im feature if the ice mantles are deposited onto the dust grains at a temperature above 
110 K (Smith et al. 1994), by as much as a factor of 3.1. And in fact, in most PPNs, the 
continuum-subtracted peak intensity of the 62 yum crystalline water-ice feature when present, 
is about a factor of ~3 lower than the 43 fim feature (see e.g. Molster ct al. 2002c). The 
continuum-subtracted peak intensity of the 43 fim feature in our spectra is 14.4 ±1.4 Jy 
(SWS 36100905) and 17.6 ± 0.8 Jy (SWS 52000931) respectively. Hence we would expect 
the ~ 62 iim feature in 1 19475 to have a peak intensity of about 4.6 Jy; given that the 
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expected width of this feature, both from laboratory and astronomical spectra, is quite large 
(about 8 /xm) this feature should have been detected in our spectra. We conclude that the 
62 nm feature is probably not present in 119475; the identification of the 43 /xm emission 
feature in our spectra as due to crystalline water ice is thus inconclusive. 

Unlike crystalline water ice, amorphous water ice produces a single broad feature with 
a peak at ~ 46 //m (Smith ct al. 1994). Molster ct al. (2002c) suggest that some amorphous 
water ice may be present in HD 161 796 in order to explain modest discrepancies between 
their model and the data in the ~40 fim wavelength region. However, in 1 19475, the 43 fim 
feature is relatively narrow, with a peak at 42.9 ± 0.04 ^m. The sharpness of this feature 
suggests that it is unlikely that amorphous water ice has a substantial presence in 119475. 
However, without detailed modelling of this feature a small contribution from amorphous 
water ice cannot be ruled out. This is because even if we assume that all ice condenses 
in crystalline form (requiring extreme conditions in the AGB outflow), UV radiation from 
the interstellar medium (ISM) can begin to amorphisize the ice (Dijkstra et al. 2003b and 
references therein). 

In the averaged SWS spectrum (Fig. Ic), a weak emission feature at ~ 23 //m attributed 
to crystalline sihcates (Molster et al. 2002b) becomes obvious. An emission feature at 48 //m 
was detected in the LWS spectrum of 1 19475 (Fig. 2). The 48 fim feature has been observed 
in several stars (see e.g. Molster et al., 2002b). Ferrarotti ct al. (2000) identified this feature 
with FeSi. However, this identification remains unconvincing since the formation of FcSi 
would require an environment where C/0 ~ 1, whereas the 48 fim feature has been observed 
in several 0-rich stars as well (Molster et al. 2002b). Furthermore, it is as yet not clear if 
this spectral feature is an artifact in the LWS instrument spectral response function, or a 
real feature (T. Lim, private communication 2006). 

The presence of [Cll] emission at 158 fim in both the on-sourcc and off-source spectra, 
combined with the absence of [01] at 63 /xm in the on-source LWS spectrum (Fig. 2) suggests 
that the [CII] emission can be attributed entirely to contamination from the ISM (see e.g. 
Castro-Carrizo et al. 2001). 
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Fig. 1. — SWS spectra of 1 19475 taken with (a) Speed 1 and (b) Speed 4 and (c) the averaged 
SWS spectrum of the star. 
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Wavelength (yum) 



Fig. 2. — (a) The observed LWS spectrum of 1 19475, (b) off-source spectrum, and (c) source 
spectrum after subtraction of the off-source data. The 48 //m emission feature seen in the 
LWS spectrum is as yet unidentified (UN) 
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Wavelength (/zm) 

Fig. 3. — The PHT-S spectrum of 1 19475 recorded with the low resolution grating spec- 
trometers PHT-SS (2.5 - 4.9 /xm) and PHT-SL (5.8 - 11.6 fim) onboard ISO. 



3.2. Spectral energy distribution (SED) 



Both SWS spectra of 1 19475 show similar flux distributions (Figs, la and b). Therefore, 
for further analysis, the two SWS observations were averaged together after rebinning the 
Speed 4 data to the resolution of the Speed 1 SWS data (Fig. Ic). The averaged SWS data 
were combined with the PHT-S data (Fig. 3) and the background subtracted LWS spectrum 
of the star (Fig. 2c). The SWS and LWS data showed good agreement in the overlap region. 
The SWS data below 14 //m and the LWS data beyond 143 //m are noisy (see Appendix A 
and Fig. 2c) and were excluded from further analysis. To reconstruct the SED, the ISO data 
were combined with the available U,B,V (Reed 2001), R,I (Monet et al. 2003), J,H,K (2Mi- 
cron All Sky Survey (2MASS); Garcia-Lario et al. 1997) magnitudes, MSX (Midcourse Space 
Experiment) fluxes and IRAS (Infrared Astronomical Satellite) fluxes of the star (Tables 2, 
3 and 4). Color-corrected IRAS fluxes (Table 4) were estimated using the IRAS Explanatory 
Supplement. We also retrieved the IRAS low resolution spectrum (LRS) from the University 
of Calgary database (http://www.iras.ucalgary.ca/~volk/getlrs_plot.html) using the "cor- 
rected raw text" option which applies the absolute calibration corrections from Cohen et al. 
(1992) to the raw IRAS LRS data. The short-wavelength (7.7-13.4 fxm) section of the LRS 
spectrum is very noisy, and has not been included. The full optical-to-far-infrared SED of 
the star is shown in Fig. 4. 




1 10 100 

Wavelength (yum) 



Fig. 4. — The observed SED of 119475. The plotted data are not corrected for interstellar 
and circumstellar reddening. U, B, V (plus symbols), R, I (triangles) fluxes are plotted along 
with J, H, K fluxes (open squares: data from Garcia-Lario et al., 1997; asterisk: 2MASS 
data), MSX (diamonds) and IRAS data (filled triangles : IRAS fluxes from the Point Source 
Catalog; filled diamonds : color-corrected IRAS fluxes). Colors highlight the PHT-S (blue), 
IRAS LRS (yellow), SWS (red) and LWS (green) data. 
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Table 2: Photometric data of 1 19475 



u 


B 


V 


R 


I 


J 


H 


K 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


10.404 


9.97 


9.312 


9.14^ 


8.78 


7.73^= 


7.41^ 


7.25^ 








9.10'' 




.- .-.TOll 
/ . / 1 o 


7. 193'' 


7.366'' 



a and b are R-band observations at epochs 1 and 2 respective^ 
2003). ^iGarcia-Lario et al (1997); 2MASS data 



y (Monet et aL, 



Table 3: MSX data 



MSX fluxes in Janskys 


Band A 
8.28 urn 


Band C 
12.13 i^m 


Band D 
14.65 iim 


Band E 
21.34 i^m 


0.214 


0.867 


2.728 


24.041 



Table 4: IRAS data 



IRAS fluxes in Janskys 


12 /im 


25 /im 


60 /im 


100 i^m 


0.54 
0.42t 


37.99 
43.9^ 


55.83 
46.7^ 


14.76 
13.5^ 



^: Color-corrected IRAS fluxes 
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4. Dust Shell Modelling of the SED 

We now describe detailed modelling of the SED of 1 19475, using the 1-D radiative 
transfer code, DUSTY (Ivezic et al. 1999). DUSTY assumes a spherical geometry and 
supports various analytical forms for the dust density distribution. The input parameters 
required for the DUSTY code are the dust temperature on the inner shell boundary (T^), 
the relative shell thickness (Y = ratio of the outer to the inner shell radius), the optical 
depth (r) at a specified wavelength, the dust grain composition and grain size distribution, 
and the spectrum of the heating radiation from the central source. 

4.1. Distance estimates 

Although the distance to the source is not a required input for DUSTY, we need it 
in order to determine the physical size and total mass of the dust shell. Various distance 
estimates for 1 19475 are available in literature. Likkel et al. (1987) used a kinematic distance 
estimate of 6 kpc. Omont et al. (1993) adopted a luminosity of 10^ Lq and estimated a 
distance of 4.8 kpc. Hrivnak & Bieging (2005) estimated a distance of 4.9 kpc to the object, 
assuming a luminosity of 8300 Lq, appropriate for a core mass of 0.63 Mq (Blocker 1995) - 
we have adopted their values for the models discussed in this paper. 

4.2. Reddening 

The total extinction towards the star, (E(B— V)=0.41), was estimated from the differ- 
ence between the observed B— V (= 0.66) value and intrinsic B— V ((B— V)o) value for the 
optical spectral type (F3Ia) of the star. For (B— V)o we used a value of +0.25 after inter- 
polating between (B— V)o for a F2Ia and a F5Ia star (Schmidt-Kaler 1982). Thus, Av=1.3, 
using Rv = 3.1, where, Rv = Av/E(B— V) (Rieke & Lebofsky 1985). Alternatively, using the 
numerical algorithm provided by Hakkila et al. (1997), which computes the 3-dimensional 
visual interstellar extinction and its error from inputs of Galactic longitude and latitude, 
and distance, from a synthesis of several published studies, we find an interstellar extinction 
of Ay — 1.5 ± 0.5. However, at the low galactic latitude (b = -1-2.73) of the star, galactic 
extinction maps (e.g. Schlegel et al. 1998) do not provide very reliable estimates of the in- 
terstellar extinction. We have corrected the observed optical and near-infrared magnitudes 
for an extinction Av=1.3, using the standard extinction laws by Rieke & Lebofsky (1985). 
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4.3. StellEir model atmosphere pcirameters 

A Kurucz model atmosphere was used as the input stellar spectrum. Several estimates 
for the effective temperature (Tcs) and gravity (log g) of 1 19475 are found in literature: Tefj 
= 7200 K, log g = 0.5 (Klochkova et al. 2002), T^g = 7750 K, log g = 1.0 (Arellano Ferro 
et al. 2001) and T^s = 7500 K, log g = 0.5 (Sivarani et al. 2001). In this paper, for the 
purpose of modeUing the SED, we adopted a Kurucz model atmosphere with Tes — 7500 K 
and log g = 1.0. 

4.4. Model Fits 

In the sections below we discuss the various model fits to the observed data. The input 
and output parameters corresponding to these models are listed in Table 5. 

The modelled flux density (reddened due to circumstellar dust) is scaled to match the 
observed flux density in the K-band (2.2 /xm) (Garcia-Lario et al. 1997), corrected for 
the total (interstellar and circumstellar) extinction. Garcia-Lario et al.'s measurement is 
consistent, within observational uncertainties, with the 2MASS K-band photometry which 
is fainter by 0.11 magnitudes. The total extinction at the K-band was estimated to be 0.14 
magnitudes. For our best fit DUSTY model (discussed later), the circumstellar extinction 
at the K-band (in optical depth units) is 0.04. 

We have chosen to use the 2.2 //m fiux rather than, e.g., an optical fiux, for scaling the 
model to the data, in order to minimise the sensitivity of our results to the following two 
modelling uncertainties. First, we do not have a very reliable estimate for the interstellar 
extinction, and second, we are using a 1-D model for an object where a significant part of the 
scattered light contribution at short wavelengths comes from non-spherical components. At 
2.2 fim, the wavelength is (i) sufficiently long that the scattered light contribution is small 
compared to the direct (but extincted) stellar contribution, and the (uncertain) extinction 
correction is small, and is (ii) sufficiently short, that the dust shell does not contribute 
significantly via emission (so our scaling is not affected by the dust properties which are to 
be derived from the model fitting). 

4.4- 1- Grain types and size distribution 

We have selected grain types for our models from the six different grain types whose 
constants are directly available in DUSTY: "warm" (Sil-Ow) and "cold" (Sil-Oc) silicates 
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from Ossenkopf et al. (1992), silicates and graphites (Sil-DL and grf-DL) from Draine 

and Lee (1984), amorphous carbon (amC-Hn) from Hanner (1988) and SiC (SiC-Pg) from 
Pegourie) (1988). We find that using sihcate grains (Sil-Ow, Sil-Oc, Sil-DL) provide good fits 
to the SED. The use of silicate grains is also supported by our discovery of silicate features 
(crystalline forsterite and pyroxenes) in the SWS spectrum of 119475 (Fig. 1). Fig. 5 shows 
the modelled fits to the observed SED using Sil-Ow, Sil-Oc and Sil-DL grains. Optical 
constants of the Sil-DL grains from Draine and Lee (1984) are based on a combination of 
laboratory measurements and astronomical observations. Optical constants of the warm 
oxygen-deficient (Sil-Ow) and cool oxygen-rich (Sil-Oc) silicates are based on observational 
determinations of the opacities of circumstellar and interstellar silicates as well as laboratory 
data (Ossenkopf et al. 1992). In estimating the optical constants of Sil-Ow and Sil-Oc, the 
effects of small spherical inclusions of various materials (oxides, sulfides, carbides, amorphous 
carbon and metallic iron) upon silicate opacities were taken into account. Fig. 5 shows that 
using Sil-Oc has the effect of producing slightly higher ffux at ~ 18 — 25;um. Sil-Ow grains 
appear to provide the best possible fit to the SED from the near to the far-infrared. We 
therefore adopted Sil-Ow for the purpose of modelling. 

The standard Mathis, Rumpl, Nordsieck (MRN) (Mathis et al. 1977) power law was 
assumed for the grain size distribution, i.e., n(a) oc a"*^ for amm <amax with q = 3.5, amin = 
0.005 /im and amax = 0.25 /xm. 
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^■^ 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

log A(/im) log X(jum) 



Fig. 5. — Model fits to the observed SED using Sil-Ow (dashed curve), Sil-Oc (dotted curve) 
and Sil-DL (dashed-dotted curve) grain types (model numbers 1, 2 and 3 in Table 5). The 
meanings of the symbols are as in Fig. 4. The PHT-S, SWS and LWS spectra are shown as 
solid curves. 
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4 ■4- 2- Density distribution 

We tried two power laws, r^^ and r^"^, for the density distribution. The best fits in the 
case of the r^^ power law for Y of 5, 10 and 15 are shown in Fig. 6. Using the r^^ power law, 
we find that even with Y as small as 5, the modelled flux systematically exceeds the LWS 
flux beyond ~70 /im. Adopting an density distribution and Y = 5, produces a good fit 
to the SED upto the limit of the LWS spectrum. For Y > 5, the modelled flux again begins 
to exceed the observed LWS flux beyond 100 //m. The best fit models using r~^ power law 
and Y = 5 and 15 are shown in Fig. 7. 
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Fig. 6. — As in Fig. 5, but showing model fits using a density distribution and relative 
shell thickness of 5 (dashed line), 10 (dashed-dotted) and 15 (dotted) (model numbers 4, 5 
and 6 in Table 5). 
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Fig. 7. — As in Fig. 5, but showing model fits using a r ^ density distribution and relative 
shell thickness of 5 (dotted) and 15 (dashed) (model numbers 7 and 8 in Table 5). 
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4.4-3. Additional constraints on the models 

Recent HST images of 1 19475 (Sahai et al. 2006) show a r'^ halo with an inner radius 
of ~ 1.5"; the surface-brightness hmited value of the outer radius is ~ 5". For our assumed 
distance (4.9 kpc), these values correspond to Rjn = 1.1 x lO^^cm and Rout = 3.6 x lO^^cm, 
which are in good agreement with our model values - e.g., in Model 7, Rin = 8.8 x lO^^cm 
and Rout = 4.4 x lO^'^'cm. 

Additional constraints on our model are provided by the observed flux of 1 19475 at 
sub- mm wavelengths. Gledhill et al. (2002) measured a flux density of 34±5.4 mJy at 850 
lim using a beamwidth of 13". The wide-band filter used by Gledhill et al. (2002) includes 
the CO J=3-2 line in its range. However, we find, using formulae from Seaquist et al. (2004) 
for calculating the fine contribution to the filter, and an estimate of the CO J=3-2 line flux 
from a multi-transition model fit to CO data on 119475, that the line emission is not likely 
to contribute more than 20-30% of the observed flux (Sahai et al. 2006). Hence the sub-mm 
flux of 1 19475 is >20 mJy. For our best-flt model (Model 7 with r^""* power law and Y=5) we 
integrated the model surface brightness distribution (in Jy arcsec^^) at 850 fim output by 
DUSTY, using a 13" diameter aperture to obtain the 850 /xm model flux density. DUSTY 
outputs the surface brightness distribution as a function of normalized (in units of Rin) radial 
distance from the central source. The radial distances were converted into angular radii (9) 
using our adopted distance of 4.9 kpc to the source. Thus Y=5 corresponds to a shell size of 
12", for which we flnd a 850 //m flux density of 10 mJy, signiflcantly lower than the observed 
value. We discuss the implications of this discrepancy for our modelling results in §5. 
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Fig. 8. — As in Fig. 5, but showing model fits with an inner shell of relative thickness 5, and 
an outer shell of relative thickness 15 (dashed) and 30 (dotted) (model numbers 9 and 10 in 
Table 5). 
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4-4-4- Dust shell mass 

The mass in the circumstellar dust shell (Mg) was computed for each of the models 
(Table 5). For an r~^ power law dust density law, 

M, = 87ri?f„(l - y-')-^Zny(Tioo/«ioo) S, (1) 

where, Rin is the inner radius of the dust shell (Table 5), Y is the relative shell thickness 
(R-out/R-in), Tioo is the shell optical depth at 100 /xm, kwo is the dust mass absorption coeffi- 
cient at 100 iim, and S is the gas-to-dust ratio. When the density distribution follows a r~^ 
power law, 

M, = 47ri?2^F(riooMoo)5. (2) 

For the models with an r~^ inner shell, and an r~^ outer shell (see § 5), 

= 47ri?2jriooMoo) S[lnYl + {Y2/Y1 - l)]/g{Yl, Y2), (3) 

where g(Yl,Y2)=(l-l/Yl2)/2 + (1-Yl/Y2)/Yl2, and Yl and Y2 are the relative shell thick- 
nesses of the inner and outer shell, respectively. The value of the dust mass absorption 
coefficient k at far- infrared wavelength is poorly known (e.g. Jura 1986). We use kiqo = 
34 cm^ gm"^, appropriate for the silicate dust used in our models (e.g. tabulated values by 
Ossenkopf at http://hera.phl.uni-koeln.de/~ossk/Jena/tables.html), and 5=200 (a typical 
value for the ejecta from AGB stars) to derive a shell mass of 1 Mq for our best-fit model 
(shell masses for this model and others are given in Table 5). The corresponding mass-loss 
rate is 3.4xlO~^M0 yr~^ {6.7x10^^ Mq yr~^) at the inner (outer) radius of the shell, assum- 
ing a constant expansion velocity of 16 km s~^ based on CO fine observations (Sahai et al. 
2006). 

Our derived shell mass should really be considered a lower limit because the model SED 
is not very sensitive to cooler dust at large radii. This is demonstrated by the fact that 
the quality of the fits to the SED from the additional models with much larger outer radii 
compared to Model 7 (fisted in Table 5), although poorer than the latter, is not drastically 
worse. 



5. Discussion 

Our model fits to the SED of 1 19475 covering the optical to the far-infrared wavelength 
range are reasonably satisfactory. The radial density distribution and physical size of the 
dust shell are also in reasonable agreement with the imaging data from HST. Since modelfing 
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of CO J=4-3 and J=2-l line emission from 119475 also suggests a r"^ radial-density shell 
(Hrivnak & Bieging 2005), the r~^ density law for I19475's dusty, molecular envelope has 
strong support. The exponent of the density law characterising the ambient circumstellar 
medium has important implications for the primary process which results in the observed 
shapes of PPNs and PNs (Sahai & Trauger 1998) - namely, the interaction of fast collimatcd 
post-AGB wind with this medium. For example, the forward shock produced by a fast wind 
expanding into an r~^ medium moves at constant velocity, but in an r~^ medium, the shock 
accelerates and is thus subject to Rayleigh-Taylor instabihties. Thus, the dynamical age 
and the morphology of the lobes in PPNs and PNs are expected to be closely related to the 
steepness of the radial density power-law distribution. 

We note however that there is a caveat to the agreement in the density law expo- 
nent between our dust model and the CO models - the latter give a mass-loss rate of 
5.2xlO~^M0 yr~^at an inner radius Rin = 3.9 x lO^^cm based on (unpublished) mid-infrared 
imaging data, which is significantly lower than our estimate at an even larger value of Rin. 
Since the CO models for inferring the mass-loss rate depend on a number of poorly-known 
quantities such as the radial kinetic temperature distribution, the inner and outer radius of 
the CO envelope, and the fractional CO abundance, we believe that our modelling of the 
dust emission provide an independent and probably more robust estimate of the shell mass 
and mass-loss rate. 

The discrepancy between the observed and model sub-mm fiux is intriguing, suggesting 
that there arc cither (1) significant quantities of cooler grains in 1 19475 than included in our 
model, and/or (2) the dust emissivity at sub-millimeter wavelengths is significantly larger 
than the value computed by DUSTY for our assumed grain type and size range. Because 
the sub-mm observations are rather insensitive to emission at radii beyond the half-power 
point of their 13" beam, extending our 12" diameter shell to larger sizes does not significantly 
enhance the model 850//m fiux, and so does not help in resolving the problem of the sub-mm 
excess. We find this expectation is supported from the results of new models in which we 
have added an r~^ envelope surrounding the r~^ density distribution in Model 7. We do 
not find any significant increase in the sub-mm wave fiux density for models with relative 
thicknesses of 15 and 30 for the outer shell (Fig. 8); moreover, in these cases the far-infrared 
modelled flux begins to exceed the observed flux. 

It is rather unlikely that the excess sub-mm fiux can be explained by the presence of 
small grains at radii inside of the half-power point of the millimeter-wave beam, which are 
sufficiently cool (i.e., for providing the excess observed submillimeter-wave fiux) because 
they are located in regions shielded from the direct stellar radiation. Such regions with 
significantly higher-than-average extinction could plausibly exist in a 3-D physical model 
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of 1 19475. However, the fundamental difficulty with such a hypothesis is that, since small 
grains are such poor emitters at sub- millimeter wavelengths, a very large mass of these grains 
is required to produce the observed submillimeter emission - and such a mass contributes 
substantially at far-infrared wavelengths (i.e., ~100 fim) for which the warm, small grains 
in our model already produce adequate emission. For example, 0.01(20 K/T(i)M0 of cold 
dust at temperature would be needed to produce the "excess" 850 iim flux of 20 mJy, 
assuming Kgso ~3 cm^ gm^^ (appropriate for small grains; Jura, Webb & Kahane 2001), 
but this dust contributes far too much emission at far-infrared wavelengths (e.g. ~20 Jy at 
100//m) to be consistent with the data. 

Although the nebula clearly has non-spherical components, simply incorporating these 
by using a 3-D model, but without the addition of large, cooler, dust grains, will not help 
resolve the problem of the sub-mm excess. We conclude that it is likely that there are 
much larger dust grains present within our model dust shell which are not accounted for by 
our adopted diist grain distribution (which follows the standard MRN prescription). Multi- 
wavelength observations of I19475's continuum flux at radio, millimeter and sub-millimeter 
wavelengths will be needed in order to probe the mass and sizes of these grains, and justify 
further quantitative modelling with modified dust grain properties. High spatial resolution 
submillimeter-wave interferometric observations with a facility like the Submillimeter Array 
could help us identify the location of these large, cool grains. One possibility is that the cool 
grains are associated with the more compact (~5") lobe-pair of this quadrupolar nebula as 
seen in the HST image (Sahai 2004, Fig. 1). 

Gledhill et al. (2002) analysed their 850/im data of 1 19475 (and other sources), together 
with the 100/x flux, assuming a fixed dust temperature of 120 K. Although their data provide 
a potentially important constraint for determining the dust content and mass in PPNs, 
their analysis and its results (dust mass, dust emissivity power-law index) are most likely 
invalid because we find from our modelling that the dust temperatures are much lower, and 
cover a large range (46-94 K). Gledhill et al.'s derived dust mass, scaled to our adopted 
distance and multiplied by our adopted gas-to-dust ratio, is O.7M0-that this number is not 
too different from our estimate of ^IMq, is purely coincidental. Our modelling underscores 
the importance of first determining the mass and temperature distribution of dust in PPNs 
carefully by fitting their full near-to-far infrared SEDs, before inferring dust masses from 
longer-wavelength data. 

There are two other wavelength regions where the modelled flux is discrepant from the 
data. First, in the ~ (7 — ll)/im range, the model flux is lower than the observed PHT-S 
data. In this wavelength region, the SED is at a minimum because both the stellar and the 
dust shell contributions are rapidly declining at wavelengths approaching this region. It is 
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likely that the observed excess flux in this region results from the presence of a small amount 
of warm dust (i.e. at temperatures larger than those in our model dust shell - e.g, 94-46 K 
in the best-fit model), and thus at radii smaller than the inner radius of the cool dust shell. 
Second, in the ~ (27 — 42)/im range, the model flux is larger than observed, suggesting 
a small excess of dust in our model shell at some intermediate temperature range within 
the full range of temperatures in our model shell. As we have discussed in Sec. 3.1, we 
believe that amorphous water-ice is not a significant contributor in this region and therefore 
it is unlikely to help in removing this discrepancy. Both discrepancies are probably best 
addressed by going to a 2-D model of the object. 



5.1. Modelling Uncertainties 

^From the off-source LWS spectrum (Fig. 2b), it is obvious that there is an increasing 
and substantial amount of background contribution to the LWS source data (Fig. 2a) beyond 
90 /im. Although considerable care was taken in subtracting the background spectrum from 
the source spectrum (Fig. 2c), this remains a possible source of uncertainty in our modelling 
- however, we do not think it affects our major conclusions. The need to adopt an r~^ (rather 
than a r~^) density law is driven by a systematic excess in the model flux above the observed 
flux over a large wavelength range, i.e. for all wavelengths longwards of ~70 fim (§4.4.2). 
This excess cannot be purely a result of the uncertainty in the background correction, since 
the total background correction upto 90 /xm is less than ~10%, with the uncertainty in the 
correction being significantly smaller. The submillimeter excess which we find in 1 19475 is 
also too large to be accounted for by uncertainties in the LWS background subtraction. 

We would like to point out that in this paper, the emphasis in the fitting procedure is not 
on the details of the mineralogy of the dust, but rather on the bulk dust component which 
is comprised of amorphous silicates. Given that we do not have a convincing detection of 
water-ice in the dust in 1 19475, a detailed investigation of models with amorphous water-ice 
(e.g., its affect on the observed 43 iim emission feature and the overall spectral shape) does 
not appear warranted. 



6. Conclusions 

We have reconstructed the spectral energy distribution of the prc-planetary nebula 
IRAS 19475-1-3119, from the optical to the far-infrared, using ISO spectra and broad-band 
photometric data. Features at 33.6 /im, 43 //m,, and possibly at 23 /im, identified in the ISO 
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SWS spectrum of the star, are most likely due to the presence of crystalline sihcates. 

The circumstellar dust envelope of the star was modelled using the 1-D radiative transfer 
code, DUSTY. We find that the dust shell has a very substantial mass, of >1Mq, in all 
models which provide a reasonable fit to the data. Our best fit is obtained for a shell with 
an r~^ radial density, inner and outer radii of 8.8x10^^ and 4.4 xlO^^cm, dust temperatures 
ranging from about 94 K to 46 K, and r(0.55 /im)=l. The mass of this shell is >1 [34 cm^ 
g^^/K(100;Lim)] [5/200]Mq, where K(100/xm) is the 100/xm dust mass absorption coefficient 
(per unit dust mass), and 6 is the gas-to-dust ratio. 

In agreement with results from optical imaging and millimeter-wave observations of CO 
emission, our model fits support an density law for 1 19475's dust shell. The exponent of 
the density law characterising the ambient circumstellar medium has important implications 
for the interaction process between the fast collimated post-AGB winds and the dense AGB 
envelopes which results in the observed shapes of PPNs and PNs. 

Our models show some discrepancies with the data. The most significant of these is that 
the observed JCMT flux at sub-millimeter wavelengths (850/im) is a factor ~2 larger than 
the model flux, suggesting the presence of dust grains in the dust shell of 119475 which are 
signiflcantly larger than those accounted for by our adopted model dust grain distribution. 

RS is thankful to NASA for financial support for this study from the Astrophysics Data 
Program (RTOP 399-20-00-08), and from the Space Telescope Science Institute through 
Program number HST-GO-09463.01. 
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Table 5: Input and output parameters for DUSTY model fits to the observed data 



Model 


Grain 


Density 


Td 




r 






Shell 


Shell Mass 


number 


type 


law 


(K) 


at 0.55 fim 


at 100 




(cm) 


thickness 


(Mo) 


1 


Sil-Ow 




95 


1.00 


9.874 X 


10-^ 


8.54 xlO^*^ 


10 


1.24 


2 


Sil-Oc 




88 


1.00 


8.404 X 


10-^ 


1.00 xlO^^ 


10 


1.45 


3 


Sil-DL 




93 


1.50 


2.631 X 


10-3 


5.70 xlO^*^ 


10 


1.47 


4 


Sil-Ow 


r-2 


96 


1.05 


1.037 X 


10-3 


8.06 xlO^^ 


5 


1.25 


5 


Sil-Ow 


r-2 


100 


1.10 


1.086 X 


10-3 


7.12 xlO^^ 


10 


2.05 


6 


Sil-Ow 


r-2 


102 


1.10 


1.086 X 


10-3 


6.71 xlO^^ 


15 


2.72 


7 


Sil-Ow 




94 


1.00 


9.874 X 


10-^ 


8.84 xl0^6 


5 


0.96 


8 


Sil-Ow 


r-3 


94 


1.00 


9.874 X 


10-^ 


8.80 xlO^^ 


15 


1.54 


9 


Sil-Ow 


T~^, r"^ 


94 


1.00 


9.874 X 


10-^ 


8.77 xlO^s 


5, 15 


2.01 


10 


Sil-Ow 


r"^, r"^ 


94 


1.00 


9.874 X 


10--^ 


8.76 xlO^^ 


5, 30 


3.63 



The inner radius of the dust shell where the dust temperature (Td) is spec- 
ified. Rin as output by DUSTY scales as L^/^ where, L is the luminosity. The 
Rin values hsted here are for a central star luminosity of 8300 Lq. 
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A. 



Data reduction 



Offline processed SWS (OLP version 10.1), LWS (OLP version 10.0) and PHT-S (OLP 
version 10.0) data were retrieved from the ISO data archive. The SWS and LWS observations 
were further processed using ISAP (ISO Spectroscopic Analysis Package) version 2.1. 



The data analysis using ISAP consisted of extensive bad data removal primarily to 
minimize the effect of cosmic ray hits. All detectors were compared to identify possible 
features. For each sub-band, 'offset' shifts were applied to bring the flux level of the 12 
detectors to a mean value (Sturm 2000). The spectra of the 12 detectors were then averaged 
using median clipping to discard points that lay more than 2.5a from the median flux. The 
averaging was done to a resolution of 300 and 1500 for SWSOl data taken with speed 1 and 
4 respectively (Table 1). The data of sub-band 3E(27.5— 29.0 //m) are generally noisy and 
unreliable (see e.g. Hrivnak et al. 2000; Hony et al. 2002). Our SWS data, below band 3A 
(16.5 /xm) is very noisy. This is also evident from the low IRAS flux at 12 fim (0.54 Jy). 
The averaged spectra of sub-band 3D (19.5—27.5 /im) were scaled in order to match the 
flux at 25 fim with the observed IRAS flux at this wavelength. Appropriate scaling factors 
were then applied to the remaining sub-bands in order to form a continuous spectrum from 
2.38-45.2 urn. 



Reduction of the LWS on and off-source observations consisted of extensive bad data 
removal using ISAP and rebinning on a fixed resolution grid of A/AA = 250. In the case of 
the on-source data, for each detector, we examined the OLP dark current estimates made 
before (Dark 0) and after (Dark 1) the observation using LWS Interactive Analysis (LIA) 
version 10.2. To form a continuous spectrum from 43—197 /xm, using ISAP_SHIFT routine, 
we applied a 'DC offset' shift to the detectors which showed different Dark and Dark 1 
values. Dark current measurements of detectors LW3 (123—152 //m), LW4 (142—171 //m) 
and LW5 (161—197 //m) showed a decreasing temporal shift. Hence a 'gain' correction was 
applied to these 3 detectors using the ISAP_SHIFT routine (Molinari) to match the flux in 
the overlapping spectral regions between the detectors. 

The background spectrum was rebinncd to the same wavelength scale as the source 
spectrum. It was assumed that the [C II] emission at 158 /im in the on-source spectrum is 



A.l. 



SWS 



A.2. 



LWS 
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entirely interstellar in origin. The background spectrum was therefore scaled in the region 
of the [C II] line (LW4 detector) to match the line strength in the on and off-source data. 
Then, using the LW4 detector as a reference the remaining detectors were scaled to form 
a continuous off-source spectrum (Fig. 2b). Subsequently, the LWS spectrum of 119475 
(Fig. 2c) was obtained by subtracting the off-source spectrum (Fig. 2b) from the on-source 
data (Fig. 2a). 



A.3. PHT-S 

The fully processed PHT-S data of 1 19475 retrieved from the ISO data archive is sci- 
entifically validated and was used directly (i.e. with no further processing). The typical 
calibration accuracy for PHT-S observations made in staring, extended source mode is bet- 
ter than ±15% (Klaas et al., 2002). The PHT-S data, lA(i) in units of surface brightness 
(W m~^ //m~^ sr~^) obtained from the archive was converted to flux density, Fx{i) units (W 
m~^ lim~^) using the receipe given in Laureijs et al. (2003) : 

F,(z) = 10-%(z)^ 

where. 



• 7;^(i) in W m ^ /im ^ sr ^ is the surface brightness for a given detector array element 

• X{i) in /im is the central wavelength of a PHT-S detector pixel; 

• ^avei^) (MJysr~^) is the average spectral response function for the conversion 
from signal in Vs~^ to surface brightness in MJysr"^; 

• ^avei^) (Vs~"^)/Jy is the average spectral response function for the conversion from 
signal in (Vs~^) to point source flux density in Jy. 

The conversion factors C^^g and C^^g (Table 6) were obtained from the PSPECAL.FITS 
file which is part of the general calibration files (Cal-G) distributed with the ISOPHOT 
Interactive Analysis (PIA), version 10.0. 
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Table 6: Conversion factors for PHT-S data 



Wavelength 

(/xm) 


(Vs-^)/(MJysr 




(Vs-i)/Jy 


Jy/(MJysr-i) 


2.46870 


3.03987 


xlO- 


-5 


0.00333517 


9.11459 


xlO- 


-3 


2.50950 


3.01950 


xlO- 


-5 


0.00328285 


9.19780 


xlO- 


-3 


2.55010 


3.23561 


xlO- 


-5 


0.00348378 


9.28764 


xlO- 


-3 


2.59070 


3.14137 


xlO" 


-5 


0.00335591 


9.36071 


xlO" 


-3 


2.63110 


3.55078 


xlO" 


-5 


0.00376118 


9.44060 


xlO" 


-3 


2.67150 


3.47183 


xlO- 


-5 


0.00364872 


9.51520 


xlO^ 


-3 


2.71180 


3.33576 


xlO" 


-5 


0.00348496 


9.57187 


xlO^ 


-3 


2.75210 


3.41233 


xlO- 


-5 


0.00354108 


9.63641 


xlO- 


-3 


2.79220 


3.62736 


xlO- 


-5 


0.00374119 


9.69574 


xlO" 


-3 


2.83230 


3.82493 


xlO- 


-5 


0.00391801 


9.76243 


xlO- 


-3 


2.87220 


4.02213 


xlO- 


-5 


0.00409953 


9.81120 


xlO- 


-3 


2.91210 


3.68381 


xlO" 


^5 


0.00373792 


9.85524 


xlO" 


-3 


2.95190 


4.81234 


xlO" 


-5 


0.00485787 


9.90628 


xlO" 


-3 


2.99170 


5.82109 


xlO" 


-5 


0.00585568 


9.94093 


xlO" 


-3 


3.03130 


7.06671 


xlO- 


-5 


0.00708768 


9.97041 


xlO" 


-3 


3.07090 


6.49100 


xlO- 


-5 


0.00648637 


1.00071 


xlO" 


-2 


3.11030 


6.87588 


xlO- 


-5 


0.00684926 


1.00389 


xlO" 


-2 


3.14970 


6.13806 


xlO- 


-5 


0.00610456 


1.00549 


xlO" 


-2 


3.18900 


7.19383 


xlO" 


-5 


0.00713846 


1.00776 


xlO" 


-2 


3.22830 


6.93928 


xlO- 


-5 


0.00687342 


1.00958 


xlO- 


-2 


3.26740 


7.03575 


xlO" 


-5 


0.00695950 


1.01096 


xlO" 


-2 


3.30650 


6.28810 


xlO" 


-5 


0.00621426 


1.01188 


xlO" 


-2 


3.34550 


6.44410 


xlO" 


-5 


0.00636565 


1.01232 


xlO" 


-2 


3.38430 


5.97321 


xlO" 


-5 


0.00590055 


1.01231 


xlO" 


-2 


3.42320 


6.31887 


xlO" 


-5 


0.00624471 


1.01188 


xlO- 


-2 


3.46190 


5.41975 


xlO- 


-5 


0.00536089 


1.01098 


xlO- 


-2 


3.50050 


5.34569 


xlO- 


-5 


0.00529456 


1.00966 


xlO- 


-2 
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Table 6: contd.... 



Wavelength 
(/iin) 


(Vs-^)/(MJysr-^) 


(Vs-i)/Jy 


Jy/(MJysr-i) 


3.53910 


4.24939x10- 


-5 


0.00421167 


1.00896 


xlO- 


-2 


3.57760 


5.09418x10- 


-5 


0.00506002 


1.00675 


xlO- 


-2 


3.61600 


4.31928x10- 


-5 


0.00430166 


1.00410 


xlO- 


-2 


3.65430 


3.78045x10" 


-5 


0.00377251 


1.00210 


xlO" 


-2 


3.69250 


2.49969x10- 


-5 


0.00250331 


9.98554 


xlO- 


-3 


3.73070 


3.35651x10- 


-5 


0.00329863 


1.01755 


xlO- 


-2 


3.76870 


3.71818x10- 


-5 


0.00364488 


1.02011 


xlO- 


-2 


3.80670 


3.82232x10- 


-5 


0.00373432 


1.02357 


xlO- 


-2 


3.84460 


3.29000x10- 


-5 


0.00320791 


1.02559 


xlO" 


-2 


3.88240 


3.63065x10- 


-5 


0.00352998 


1.02852 


xlO- 


-2 


3.92010 


3.44827x10- 


-5 


0.00334782 


1.03000 


xlO- 


-2 


3.95780 


3.76135x10- 


-5 


0.00364328 


1.03241 


xlO- 


-2 


3.99540 


3.25984x10- 


-5 


0.00315459 


1.03336 


xlO- 


-2 


4.03280 


3.26180x10- 


-5 


0.00315082 


1.03522 


xlO- 


-2 


4.07020 


2.89407x10- 


-5 


0.00279106 


1.03691 


xlO- 


-2 


4.10760 


3.27721x10- 


-5 


0.00315620 


1.03834 


xlO- 


-2 


4.14480 


2.82243x10- 


-5 


0.00271870 


1.03815 


xlO- 


-2 


4.18190 


3.10457x10- 


-5 


0.00298790 


1.03905 


xlO" 


-2 


4.21900 


2.72054x10" 


-5 


0.00261666 


1.03970 


xlO" 


-2 


4.25600 


2.69869x10- 


-5 


0.00259448 


1.04017 


xlO- 


-2 


4.29290 


2.27141x10" 


-5 


0.00218314 


1.04043 


xlO- 


-2 


4.32970 


2.39768x10" 


-5 


0.00230451 


1.04043 


xlO- 


-2 


4.36640 


2.07467x10- 


-5 


0.00199824 


1.03825 


xlO- 


-2 


4.40310 


2.35388x10- 


-5 


0.00226821 


1.03777 


xlO- 


-2 


4.43970 


1.94929x10" 


-5 


0.00187972 


1.03701 


xlO- 


-2 


4.47610 


2.18277x10" 


-5 


0.00210694 


1.03599 


xlO- 


-2 


4.51250 


1.81484x10" 


-5 


0.00175392 


1.03473 


xlO- 


-2 
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Table 6: contd.... 



Wavelength 
(/iin) 


(Vs-^)/(MJysr-^) 


(Vs-i)/Jy 


Jy/(MJysr-i) 


4.54890 


1.98181x10- 


-5 


0.00191809 


1 


03322 


xlO" 


-2 


4.58510 


1.62796x10- 


-5 


0.00157838 


1 


03141 


xlO" 


-2 


4.62120 


1.42384x10- 


-5 


0.00138326 


1 


02934 


xlO- 


-2 


4.65730 


1.44786x10" 


-5 


0.00140976 


1 


02703 


xlO- 


-2 


4.69330 


1.50828x10" 


-5 


0.00147225 


1 


02447 


xlO" 


-2 


4.72920 


1.36050x10" 


-5 


0.00132851 


1 


02408 


xlO" 


-2 


4.76500 


1.80055x10" 


-5 


0.00176359 


1 


02096 


xlO" 


-2 


4.80080 


1.34228x10" 


-5 


0.00131913 


1 


01755 


xlO" 


-2 


4.83640 


1.44308x10" 


-5 


0.00142338 


1 


01384 


xlO" 


-2 


5.83960 


12.8575x10- 


-5 


0.0115198 


1 


11612 


xlO- 


-2 


5.93380 


11.4284x10- 


-5 


0.0102716 


1 


11262 


xlO" 


-2 


6.02800 


11.2309x10- 


-5 


0.0100773 


1 


11448 


xlO" 


-2 


6.12200 


11.6475x10- 


-5 


0.0104899 


1 


11035 


xlO" 


-2 


6.21600 


10.9303x10- 


-5 


0.00983196 


1 


11171 


xlO" 


-2 


6.30990 


8.64221x10- 


-5 


0.00776619 


1 


11280 


xlO" 


-2 


6.40380 


10.6843x10- 


-5 


0.00964404 


1 


10787 


xlO" 


-2 


6.49750 


9.88445x10- 


-5 


0.00891797 


1 


10837 


xlO" 


-2 


6.59120 


9.40491x10" 


-5 


0.00848214 


1 


10879 


XlO" 


-2 


6.68480 


8.12382x10" 


-5 


0.00732521 


1 


10902 


XlO" 


-2 


6.77830 


9.59879x10" 


-5 


0.00870225 


1 


10302 


XlO" 


-2 


6.87170 


9.17306x10- 


-5 


0.00831815 


1 


10278 


XlO" 


-2 


6.96500 


10.4062x10" 


-5 


0.00944026 


1 


10232 


xlO" 


-2 


7.05830 


9.92587x10" 


-5 


0.00900942 


1 


10172 


xlO" 


-2 


7.15150 


9.66050x10" 


-5 


0.00877220 


1 


10126 


xlO" 


-2 


7.24460 


8.08621x10" 


-5 


0.00734715 


1 


10059 


xlO" 


-2 


7.33760 


7.66112x10" 


-5 


0.00696665 


1 


09968 


xlO" 


-2 


7.43060 


6.90864x10" 


-5 


0.00628884 


1 


09856 


XlO" 


-2 



-33- 



Table 6: contd.... 



Wavelength 
(/iin) 


(Vs-^)/(MJysr-^) 


(Vs-i)/Jy 


Jy/(MJysr-i) 


7.52350 


8.15585x10- 


-5 


0.00743328 


1 


09721 


xlO- 


-2 


7.61630 


7.35407x10- 


-5 


0.00671185 


1 


09568 


xlO" 


-2 


7.70900 


7.07851x10- 


-5 


0.00647062 


1 


09395 


xlO- 


-2 


7.80160 


6.83308x10" 


-5 


0.00625740 


1 


09200 


xlO" 


-2 


7.89420 


7.86450x10" 


-5 


0.00716019 


1 


09836 


xlO" 


-2 


7.98660 


7.69969x10" 


-5 


0.00702496 


1 


09605 


xlO" 


-2 


8.07900 


7.75643x10" 


-5 


0.00709291 


1 


09355 


xlO" 


-2 


8.17140 


6.57546x10" 


-5 


0.00602763 


1 


09089 


xlO" 


-2 


8.26360 


6.08753x10" 


-5 


0.00559467 


1 


08809 


xlO" 


-2 


8.35580 


5.63152x10- 


-5 


0.00514833 


1 


09385 


xlO- 


-2 


8.44780 


6.16770x10- 


-5 


0.00565487 


1 


09069 


xlO- 


-2 


8.53980 


5.88137x10" 


-5 


0.00540839 


1 


08745 


xlO- 


-2 


8.63180 


6.33101x10" 


-5 


0.00583982 


1 


08411 


xlO" 


-2 


8.72360 


3.26216x10" 


-5 


0.00299513 


1 


08915 


xlO" 


-2 


8.81540 


5.01120x10" 


-5 


0.00490196 


1 


02228 


xlO" 


-2 


8.90700 


5.79664x10" 


-5 


0.00562173 


1 


03111 


xlO" 


-2 


8.99860 


5.44365x10" 


-5 


0.00519205 


1 


04846 


xlO" 


-2 


9.09020 


5.10184x10- 


-5 


0.00483174 


1 


05590 


xlO- 


-2 


9.18160 


4.59752x10- 


-5 


0.00428244 


1 


07357 


xlO- 


-2 


9.27300 


4.28071x10" 


-5 


0.00396460 


1 


07973 


xlO" 


-2 


9.36430 


4.05319x10" 


-5 


0.00373380 


1 


08554 


XlO" 


-2 


9.45550 


4.54677x10" 


^5 


0.00412516 


1 


10229 


xlO" 


-2 


9.54660 


4.58967x10" 


-5 


0.00414712 


1 


10671 


xlO" 


-2 


9.63760 


3.91220x10" 


-5 


0.00348406 


1 


12289 


xlO" 


-2 


9.72860 


3.55302x10" 


-5 


0.00315427 


1 


12642 


xlO" 


-2 


9.81950 


2.83132x10" 


-5 


0.00247981 


1 


14175 


xlO" 


-2 


9.91030 


3.39671x10- 


-5 


0.00296851 


1 


14425 


xlO" 


-2 
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Table 6: contd.... 



Wavelength 
(//m) 


(Vs-^)/(MJysi 




(Vs-i)/Jy 


Jy/(MJysr-^ 


) 


10.0010 


3.34585x10- 


-5 


0.00288679 


1 


15902 


xlO" 


-2 


10.0917 


2.71791x10- 


-5 


0.00234248 


1 


16027 


xlO- 


-2 


10.1823 


2.93380x10- 


-5 


0.00249772 


1 


17459 


xlO- 


-2 


10.2727 


3.21568x10- 


-5 


0.00273768 


1 


17460 


xlO- 


-2 


10.3632 


2.83052x10" 


-5 


0.00238194 


1 


18833 


xlO- 


-2 


10.4535 


2.61391x10- 


-5 


0.00220172 


1 


18721 


xlO- 


-2 


10.5437 


2.52504x10- 


-5 


0.00210409 


1 


20006 


xlO" 


-2 


10.6339 


2.60993x10" 


-5 


0.00217871 


1 


19792 


xlO" 


-2 


10.7240 


2.67050x10" 


-5 


0.00220733 


1 


20983 


xlO" 


-2 


10.8140 


2.50139x10" 


-5 


0.00207279 


1 


20677 


xlO" 


-2 


10.9040 


2.50909x10" 


-5 


0.00206055 


1 


21768 


xlO" 


-2 


10.9938 


2.42555x10" 


-5 


0.00199853 


1 


21367 


xlO" 


-2 


11.0836 


2.15947x10" 


-5 


0.00176475 


1 


22367 


xlO" 


-2 


11.1733 


2.05535x10- 


-5 


0.00166620 


1 


23356 


xlO" 


-2 


11.2629 


1.76285x10- 


-5 


0.00143567 


1 


22789 


xlO" 


-2 


11.3524 


1.74882x10- 


-5 


0.00141409 


1 


23671 


xlO" 


-2 


11.4419 


2.05966x10" 


-5 


0.00167428 


1 


23018 


xlO" 


-2 


11.5313 


2.31557x10" 


'5 


0.00187026 


1 


23810 


xlO" 


-2 


11.6206 


1.98068x10" 


-5 


0.00159024 


1 


24552 


xlO" 


-2 
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